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SUMMARY 


The spectral absorption properties of C 3 have been measured in a shock 
tube containing a test-gas mixture of acetylene diluted with argon. The 
absorption of a pulsed xenon light source was measured by means of eight photo- 
multiplier channels to a spectrograph and an accompanying drum camera. The 
postshock test-gas temperature and pressure were varied over the ranges 3240 to 
4300 K and 37 to 229 kPa, respectively. 

The results showed appreciable absorption by C 3 for the wavelength range 
300 to 540 nm. The various reported measurements of the heat of formation of 
C 3 which are available in the open literature have been reviewed, and a value 
of 198 kcal/mol is recommended. This value, along with best available values 
for other species, was used to calculate the number density of C 3 for the con- 
ditions of the present experiments in order to compute absorption cross section 
or electronic oscillator strength. The computed electronic oscillator strength 
varied from a high of 0.062 at 3300 K to a low of O.O 36 at 3900 K, 


INTRODUCTION 

Entry into the atmospheres of the outer planets is accompanied by large 
radiative heating of the entry probe forebody. To some degree, this thermal 
radiation is blocked (absorbed) by the cooler gases being evolved from the 
ablative heat shield. For example, reference 1 shows that the ablation gases 
emitted from a phenolic-carbon heat shield absorb about half the incident radi- 
ant energy for an entry into the atmosphere of the planet Jupiter. Although 
the calculation method used in reference 1 is the best available estimate of 
the absorption effects of the ablator gases, its accuracy is hampered by the 
limited knowledge of the absorption properties of the various gas species pres- 
ent, Under the assumption of thermochemical equilibrium, the calculations pre- 
dict that a significant fraction of the ablator gas will be composed of such 
molecules as C 3 , C 4 , C 5 , C 2 H, and C 3 H, for which little or no information is 
available on spectral absorption properties. 

Reference 1 indicates that C 3 is the dominant species near the ablator 
wall, with a mole fraction exceeding 0,5 for some cases. Thus, an accurate 
determination of its spectral absorption cross section will greatly benefit 
the radiant transfer calculations. The present study examines the absorption 
by triatomic carbon C 3 in the spectral range 260 to 560 nm. 

The most prominent spectral feature of C 3 is a continuum-^and system in 
the vicinity of 400 nm. Frequently referred to as the "X4050 A bands" 

(ref. 2 ), and first observed in cometary spectra, this line group has been 
shown (ref. 3 ) to be due to a ground-state-connected electronic transition in 
the C 3 molecule. Only a few lines are observed at cryogenic temperatures 
(e.g., refs. 4 to 7), but at elevated temperatures, a strong apparent contin- 
uum dominates this spectral region (ref. 2). This apparent continuum was shown 


by Brewer and Engelke (ref. 8 ) to be due to the same bound states as the 
observed lines and, therefore, not a true continuum but a "pseudocontinuum” of 
very closely spaced rotational lines of the electronic transition. Reference 8 
includes an absorption measurement in saturated carbon vapor at a temperature 
of 3200 K. The estimate of C 3 number density given in reference 8 was used in 
reference 9 to convert these data to spectral cross section. These appear to 
be the only quantitative data on C 3 absorption at high temperatures that are 
available in the open literature. 

In the present study, the species C 3 was produced by shock heating an 
acetylene-argon test-gas mixture in a shock tube. The calculated temperature 
behind the incident shock wave, based on shock velocity measurements, ranged 
from about 3240 to 4300 K. The postshock pressure varied from 37 to 229 kPa. 

The postshock test gas was viewed in absorption by means of a high-pressure 
xenon pulsed light source. The number density was calculated by using the mea- 
sured shock velocity and assuming thermochemical equilibrium. From the measured 
absorption and the calculated number density, the spectral absorption cross sec- 
tion was calculated for the spectral range 260 to 560 nm. The variation of this 
spectral absorption cross section with temperature is obtained for the tempera- 
ture range of the experiments . 


SYMBOLS 


Values are given in SI Units except for heat content, which is given in 
calories (1 cal = 4.184 J) . Following common practice, the abbreviation eu is 
used for entropy units of calories per degree kelvin per mole. 


a 


30*^1 

a2,a3 


c 


Cp 

D 

E 

e 

e(t) 


F(t) 

F’(t) 


f 


radius of carbon particulates, nm 

coefficients in curve fit of absorption cross section as a function 
of temperature (eq. ( 8 )) 

velocity of light, cm/sec 

specific heat at constant pressure, cal/mol-K 

film density 

film exposure, arbitrary units 

charge of an electron, C 

time-dependent emission of the test gas as measured by a 
spectrograph-channel photomultiplier tube 

time-dependent absorption defined by equation ( 1 ) 

absorption parameter corrected for gas emission, as defined by 
equation ( 2 ) 

focal length, m 
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fe electronic oscillator strength 

Gx Gibbs free energy at temperature T, cal/mol 

h static enthalpy, cal/mol 

I light intensity, arbitrary units 

Iq unattenuated light intensity, arbitrary units 

Ig light intensity after passing through test gas with path length s, 

arbitrary units 

extinction coefficient at wavelength X or wave number o), cm"^ 
m mass of electron, kg 

in mass of carbon particulates per unit flow volume, kg/m^ 

N number density, cm“3 

N-|,N 2 real and imaginary components, respectively, of refractive 
index of carbon 

p pressure , Pa 

P;^ Rayleigh factor 

R universal gas constant, 8314 J/kmol-K (1.98? cal/mol-K) 

Sx entropy at temperature T, cal/mol-K 

S(t) time-dependent photomultiplier signal 

s optical path length through test gas, cm 

T temperature , K 

T 2 calculated equilibrium test-gas temperature immediately behind 

incident shock wave, K 

t time , sec 

Vg velocity of incident shock wave, km/sec 

X distance along shock tube measured from diaphragm, m 

Y slope of linear portion of film sensitivity curve (eq. (3)) 

Y mean isentropic exponent for test gas 
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AHf 


heat of formation at pressure of 1 atm (101 kPa) and temperature T, 
kcal/mol 


X 

Pc 

c^X ’*^0) 


wavelength, nm 
density of graphite, kg/m3 

absorption cross section at wavelength X or wave number w, cm^ 
0 ) wave number, om“^ 

Subscripts : 
abs absorption 

p particles 

sea scattering 

0 prerun condition for light signal 

1 condition before shock arrival 

2 condition behind incident shock 

50 calculated conditions 50 ysec after shock passage 


EXPERIMENTAL APPARATUS 
Shock Tube 

The shock tube used in the present investigation is the same facility 
described in reference 10. It is an all-stainless-steel tube 15.2 cm in 
internal diameter with a driven tube approximately 15.4 m long. The electric 
arc electrode assembly was removed from the driver chamber and a driver exten- 
sion section was added to make a driver chamber length of 2.44 m. The boundary- 
layer splitter-plate assembly in the test section was not used, since in the 
present study the postshock flow Mach numbers were less than 2, and preliminary 
tests indicated that the splitter-plate assembly sometimes choked the flow. 

Thus, the test-section windows were mounted flush with the tube walls. 

Helium was used as the driver gas for all runs. After preliminary evacua- 
tion, helium was injected into the driver chamber until the diaphragm ruptured. 
Diaphragms were stainless steel with the conventional grooves machined in a 
cross pattern to cause the diaphragm to open in the usual four-petal 
configuration. 

The driven chamber was initially evacuated to a pressure of approximately 
0.03 Pa by using a diffusion pump and a liquid-nitrogen-cooled trap. The test 
gas was injected into the chamber approximately 3 to 5 minutes prior to the run. 
The increase in chamber pressure due to outgassing and leakage was approximately 
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0.07 Pa/min. Since the initial test-gas pressure was 700 Pa or greater, the 
ratio of air and water vapor to test gas was believed to be 6 x 10”^ or less 
for all runs . 


Optical Instrumentation 

A schematic drawing of the optical arrangement and shock tube is shown in 
figure 1 . The light source was a high-pressure xenon lamp which was pulsed by 
the discharge of a capacitor-inductor transmission line to produce a nearly 
constant light intensity for about 1 msec. This pulse was initiated by the 
passage of the shock wave at a station upstream of the test section, so that 
the light pulse was initiated 200 to 300 psec prior to the arrival of the shock 
wave at the window location. Quartz lenses were used to form an image of the 
lamp on a small orifice 1.8 mm in diameter and then to form a collimated beam 
using the image at the orifice as the apparent source. Another lens focused 
this beam on the spectrograph entrance slit. Field stops 1 cm in diameter were 
used at the window locations and along the collimated beam path to reduce stray 
light and define the geometry of the optical path. 

The spectrograph used was a Czerny-Turner mount with a nominal aperture of 
f/6.3. Internal masks were used to insure that the spectrograph was filled 
with light. These masks limited the spectrograph to approximately f/8. The 
spectrograph was modified by inserting a front-surface mirror which reflected 
approximately 80 percent of the light into a rotating-mirror camera which had 
been adapted to the spectrograph. The remaining 20 percent of the light passed 
over the top of the mirror and fell on the focal plane. Eight light pipes were 
mounted in a group with their entrance faces at the focal plane. Each light 
pipe led to a photomultiplier tube. When used in the ultraviolet range, the 
entrance faces of the light pipes were coated with sodium salicyclate. The 
effective width of the light pipes ranged from 1.1 to 1.9 mm. Therefore, the 
effective band pass was 4.3 to 7.7 nm and 2.1 to 3-8 nm for the two gratings 
used in the investigations. The gratings used were a 300 line/mm grating 
blazed at 500 nm and a 600 line /mm grating blazed at 300 nm. The approximate 
values of linear dispersion at the focal plane were 41.0 nm/cm for the 
300 line/mm grating and 20.2 nm/cm for the 600 line/mm grating. 

The photomultipliers were used with load resistors that limited output 
current to less than 1 mA. The linearity of each photomultiplier tube was 
checked with neutral-density filters. An example of the linearity check is 
shown in figure 2 for one channel. On each channel, ballast capacitors paral- 
leled each of the last two resistors of the dynode chain so as to maintain 
nearly constant dynode voltages for extended signal periods. 

The drum camera was a commercially manufactured device which was modified 
to adapt to the spectrograph. It used an air-driven turbine to turn the drum 
318 ± 5 rps. The length of film (internal circumference of the drum) was 
31.88 cm, and therefore the film speed was 10.1 ± 0.1 cm/msec. Since the 
effective slit height at the film was approximately 1.5 mm, the time resolu- 
tion on the film was no better than 15 psec. 
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Velocity-Pressure Instrumentation 

The distances from the primary diaphragm to the centers of the instrumen- 
tation ports are shown in table I. A thin- film resistance thermometer was 
installed in station 1 to detect the passage of the shock wave. The output 
signal from this gage was used to initiate both the system of time-interval 
meters which measured the shock transit times and the time-delay generators 
which triggered the oscilloscopes. 

The remaining instrumentation ports had either thin-film gages, light 
pipes leading to photomultipliers, or quartz pressure transducers mounted as 
indicated in table I. The thin-film gages were used because they gave low- 
noise, rapid-rising signals which were easily amplified for use as trigger sig- 
nals. The photomultipliers merely gave an indication of the approximate test 
time. The pressure transducers showed the steadiness of the test-gas slug. 
Examples of the oscillograms from a photomultiplier and a pressure transducer 
are shown in figure 3 . 

The shock velocity was determined from the measurements of time-interval 
meters. The amplified signal from the instrumentation mounted at each station 
was used to stop a time-interval meter. The start gate for all the time- 
interval meters was provided by the thin-film gage at station 1 . An example 
of the calculated average velocity between successive stations is shown in fig- 

m/sec 

ure 4. The shock attenuation ranged from 5 to 20 depending on the pres- 

m 

sure level. The velocity which was used to reduce the data was the final veloc- 
ity at the test section. 


Test Gas 

The test gas was purchased already mixed and certified by the supplier to 
be 2.99 percent C 2 H 2 , 97.01 percent argon by volume. A few runs were made as 
a check with a 1.97 -percent C 2 H 2 , 98.03-percent argon mixture. Twice during 
the test program, a sample of the test gas was analyzed by gas chromatography 
to insure that the mixture remained uniform. 


Test Procedure 

It was found that if the lamp had not been fired for several hours, it was 
necessary to fire it several times in order to insure reproducible light pulses. 
After these warmup firings, the lamp regularly reproduced an identical pulse 
within ±5 percent. When the lamp had been fired sufficiently to Insure good 
repeatability, a prerun shot was recorded from the eight photomultipliers con- 
nected to the spectrograph light pipes. These records were regarded as the 
unabsorbed intensity as a function of time for each channel. The capacitor 
bank which fired the lamp was charged to 355 V, as read on a digital voltmeter, 
for the prerun record. 
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As mentioned previously, the driven chamber was evacuated to about 0.03 Pa 
prior to each run. The driver chamber was evacuated to about 1.5 Pa by using a 
mechanical pump only. After the test gas was injected into the driven tube, 
the lamp supply capacitors were charged to approximately 365 V. Since the 
charge slowly leaked off the capacitors, the succeeding run operation steps 
were paced in an attempt to synchronize the diaphragm rupture with a capacitor 
voltmeter reading of 355 V. This result was accomplished within ±2 V for 
nearly all runs. 

Following each run, the shock-tube windows were cleaned, since postrun 
lamp firings showed that the intensity was reduced by 10 percent or so as a 
result of deposits on the windows. The entire shock-tube interior was cleaned 
periodically (approximately every 20 runs) with ethyl alchohol. 


DATA REDUCTION PROCEDURES 
Spectrograph Channels 


A sample record from a spectrograph channel is shown in figure 5. Fig- 
ure 5(a) shows the prerun record of the light pulse, and figure 5(b) shows the 
oscillogram taken during the run. The arrival of the shock wave is clearly 
indicated by the sudden drop in photomultiplier signal level. Figure 5(c) 
shows the record after being digitized, along with the digitized prerun record 
for the same channel. Figure 5(d) presents the fraction of light absorbed 


F(t) 


So(t) - S(t) 
So(t) 


( 1 ) 


where S(t) and Sg(t) are the time-varying values of the photomultiplier 
signal for the run and prerun, respectively. For these traces, time equal to 
zero corresponds to the start of the light pulse. 

The parameter F(t) was plotted for each channel of each run by a program- 
mable calculator with the digitized oscillograms as input data. Typically, the 
value of F(t) increased moderately with time after shock passage, most likely 
because of the growing wall boundary layer. Many runs showed that the parameter 
F(t) took 40 psec or so to reach a plateau value, suggesting that the chemical 
relaxation rate took approximately this long to adjust to the postshock 
condition. 


The light received by the spectrograph included the emission from the gas 
as well as the (partially absorbed) emission from the lamp. In order to deter- 
mine the contribution from gas emission, runs were made with the lamp turned 
off. For some runs, a light chopper was employed so that the absorption and 
emission could be determined on the same run. The field stops in the light 
path attempted to insure that the volume of gas contributing to the emission 
was essentially the same as that for which the absorption was determined. 


The emitted light was expressed as a fraction of the prerun lamp signal 
Sq, and this quantity was used to correct the absorption parameter F(t). If 
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the corrected absorption is denoted as F'(t) and the emission as e(t), the 
following expression results: 


F’(t) 


So(t) - [S(t) - e(t)] e(t) 

= F(t) + 

So(t) So(t) 


( 2 ) 


For most run conditions, the emission term was immeasurably small. For the 

e(t) 

highest temperature runs, the emission term grew as large as ^ =0.1 and 

represented a correction to the value of F(t) of the order of 30 percent. A 
few runs were made in pure argon test gas to insure that the emission and 
absorption from the carrier gas were negligible. 


Drum Camera Film 

A sample of the film obtained from the drum camera is shown in figure 6(a). 
The prominent features of the absorption are an apparent continuum and a few 
bands. All the bands that were observed over the entire wavelength range of the 
study belong either to C2 or the violet system of CN. This latter system is so 
strong that only a small contamination of air in the test-gas mixture causes it 
to appear. 

Figure 6(b) shows a microdensitometer tracing of the film for the same run 
as that of figure 6(a). Two traces are shown, one taken at a time just before 
shock arrival and the other taken about 100 ysec after shock arrival. The dif- 
ference between the traces is a measure of the absorption. 

The film yields so much more spectral detail than the light pipe channels 
that an attempt was made to use the densitometer traces quantitatively. The 
following procedure discusses how this was done. It must be emphasized, how- 
ever, that only limited confidence can be attached to the quantitative film 
density results, because of the many pitfalls associated with using film quan- 
titatively and the approximate procedure used in the present study. 

Twice during the investigation, sensitometer samples were made on film 
samples and these were developed in the same manner as the test films. The 
sensitometer samples were then run on the microdensitometer and the specular 
density of the film was plotted against the logarithm of the exposure, as shown 
in figure 7. The linear portion of this curve showed a slope of approximately 
Y = 0.9. The film density D was thus assumed to be related to exposure by 


D-i - D 2 = y[iogio (E-|) - log-10 (^2)] 


(3) 


where the subscripts 1 and 2 indicate times just before and just after 
shock arrival, respectively. Exposure is defined as the product of the light 
intensity and time; therefore, equation (3) may be rewritten 
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10(Di-I>2)A 


(4) 


E-i I-|t 

E2 l2fc 

Since the exposure time is the same for each element on the test film, 
equation (4) yields the ratio of light intensity before shock arrival to that 
after the shock as a function of density difference on film. This ratio was 
then used to determine the wavelength-dependent extinction coefficient K^. As 
noted previously, the values of extinction coefficient determined from the film 
must be viewed with great caution. However, as is shown subsequently, they show 
quite good repeatability and are in essential agreement with the data obtained 
from the light pipe-photomultiplier channels. The two most daring assumptions 
in the procedure just described appear to be 

( 1 ) The determination of the slope y was made by using the xenon flash 
lamp in the sensitometer , and these "white” light data were then assumed to 
apply for all wavelengths throughout the spectral range of the investigation 
and for all film exposures, even though at low film density, the slope was 
distinctly lower. 

(2) The same calibration data were assumed to apply for all runs; thus any 
differences in film development from run to run. were ignored. 


Absorption Cross Section 

The extinction coefficient K\ is defined by the relation 

dl = -Kx(s)I ds (5) 

(See ref. 11 for a review of fundamentals in absorbing media.) If the absorbing 
medium has uniform absorbing and scattering properties along the length of the 
light path s, equation (5) may be integrated from 0 to s to yield 

lo 

In — 

Is 

Kx = (6) 

s 

Here, Ig and Ig denote, respectively, the intensity as the light beam enters 
and leaves the absorbing medium of length s. 

In the present study, scattering effects were ignored and the absorption 
coefficient was assumed identical with the extinction coefficient. Solid parti- 
cles, if present, would be expected to be a significant scattering agent. For 
this reason, the present test conditions avoid the lower temperature range, 
where solid carbon particles might be expected to be present. In "Results and 
Discussion,” tests are described which support the conclusion that the present 
experiments are not significantly affected by solid particle scattering. 
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An absorption cross section may be defined as 


KX 



(7) 


Use of the absorption cross section to define the absorption properties of a 
gas accounts for the dependence of absorption coefficient on gas density by 
assuming a linear variation with the number of molecules available to absorb. 
Reduction and use of the present data therefore require that the number density 
of Cg molecules be determined. 

For the present tests , the number density of molecules was determined 
by a free-energy-minimization calculation assuming thermochemical equilibrium 
in the test region behind the shock wave. The chemical equilibrium computer 
program ACE of reference 12 was used for this calculation; the thermodynamic 
state properties ahead of the shock wave and the shock velocity for each run 
were required inputs to the program. The program also requires a value for 
heat of formation and a definition of the temperature-dependent thermodynamic 
variables for each gas species which is to be considered in the final mixture. 
This requirement creates a problem in that the heat of formation of has not 
been determined as accurately as is desirable. Measurements of the heat of 
formation of range from AHq = 188 kcal/mol (ref. 13) to as high as 
AH298 = 212 kcal/mol (ref. 14). This spread in heat of formation, together 
with the uncertainty in thermal properties, produced an unacceptably large 
difference in the calculated number density of C^. Consequently, a review of 
the various measurements of the heat of formation and thermal property calcula- 
tions of C3 was included in the present study and is reported in the appendix. 
The thermal property calculations of Lee and Sanborn (ref. 15) and a heat of 
formation AHq = 198 kcal/mol were selected and used in the ACE program for 
computing the C^ number density. 


Varying Stream Properties 

As was stated previously, the absorption parameter F'(t) usually 
increased moderately with time after shock passage. Three possible causes of 
this time variation of absorption are suggested: (1) the flow might not be in 

chemical equilibrium; (2) the free-stream properties vary as a result of the 
displacement effects of the boundary layer; and (3) the boundary layer, which 
grows with distance behind the shock wave, absorbs more per unit optical path 
length than does the inviscid flow. In this section, these possible causes are 
examined in turn in an effort to evaluate their effects on the present data. 

Thermochemical equilibrium .- Since rate constants for the important reac- 
tions are not well known, a finite-rate calculation of the reaction process 
would not be reliable. However, several facets of the present data suggest that 
chemical equilibrium was essentially obtained in about 40 psec after shock pas- 
sage. (It might be noted that the residence time of the particles was about 
4 times greater than the observation time. Thus, at an observation time of 
40 ysec, the observed gas particles passed through the shock wave approximately 
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160 ysec earlier.) These various facets that suggest chemical equilibrium for 
the present test conditions are now briefly enumerated: 

(1) The photomultiplier tubes used to detect shock wave and interface pas- 
sage typically displayed signals which rose to a level in 20 to 40 psec. There- 
after, the signal level either decayed somewhat (as in the example in fig. 3(a)) 
or remained essentially constant. Since these light signals were not filtered 
spectrally, they contained contributions from all emitters within the spectral 
range of the detector (possibly including impurities) and would be expected to 
be quite sensitive to temperature variations. Thus, the relative flatness of 
the signal would seem to imply only small temperature variation with time. 

(2) On runs in which the emission was measured with the spectrograph (i.e., 
the absorption lamp source was not used), the emission at each wavelength rose 
to a plateau in about 20 Usee after shock passage and then remained remarkably 
constant for the remainder of the test time. This emission was assumed to be 
due to C3 and therefore implied constant number density and temperature of the 
C3 species. 

(3) Timewise surveys of the drum camera films at the wavelengths of the 
band heads of the C2 Swan band system were also very constant with time. Thus, 
the absorption due to C2 was essentially time invariant. Also, the measured 
absorption coefficient was in fairly good agreement with that predicted by 
reference 16. 

(4) The results of runs for a test-gas mixture of 2 percent acetylene in 
argon were in good agreement with the results for a 3-percent acetylene mixture. 
Also, no significant trends with pressure level were observed. 

Effects of varying free-stream properties .- Mirels has shown (ref. 17) 
that the free-stream flow properties vary with distance behind the shock wave 
because the growing wall boundary layer changes the effective cross-sectional 
area of the inviscid flow. Reference 17 presents relationships for calculating 
the varying flow properties, for either wholly laminar or wholly turbulent 
boundary-layer growth, based on the assumption of an isentropic process in the 
postshock inviscid flow. For the present data, a turbulent boundary layer has 
been assumed because of the Reynolds number range of the test conditions. This 
assumption appears justified, since the measured test times were in good agree- 
ment with those calculated by the method of reference I8 for a turbulent bound- 
ary layer. 

Four runs were selected as approximately bounding the pressure-temperature 
range of the present tests . An isentropic exponent based on the equilibrium 
mole fractions of the major species was calculated for each of these runs. The 
mean value of this isentropic exponent y was 1.654, and this value was used 
for all runs to compute the isentropic compression of the shocked gas, starting 
with the postshock equilibrium conditions. Figure 8 shows the computed temper- 
ature variation at the test location for the four sample runs. It may be seen 
that the temperature variation T/T2 is very similar for all four runs; there- 
fore, for all runs of the test program, the time variation of the postshock 
temperature was modelled as 
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T/T2 = 1 + 0.0066(t X 1o6)0.3 

This relation is also shown in figure 8. The pressure variation was then found 

Y / (y - 1 ) 

by the isentropic relation p <x T*'^ ' 

In a similar manner, the variation of the mole fraction of C3 was curve 
fit for the present test range as a function of pressure and temperature. This 
permitted determination of the time-varying number density as a perturbation of 
the value calculated immediately behind the shock wave. 

Figures 9(a), 9(b), and 9(c) show examples of the effect of applying the 
time-varying number density to the absorption data . As illustrated by figure 9 , 
this calculation method was not uniformly successful in accounting for the time 
variation of the absorption. In fact, for the highest temperature conditions, 
the absorption cross section varies more with time than does the absorption 
factor F'(t). It seems apparent that correction for the variation of free- 
stream properties will not by itself yield a time-invariant absorption cross 
section. 

Boundary-layer absorption .- No analysis of the absorption within the 
boundary layer was attempted in the present study, and for the data reduction 
a constant value for the absorption path length equal to the tube internal 
diameter was used. This would imply that the number density of C3 per unit 
path length in the boundary layer was the same as that in the inviscid stream. 
Such a case would not usually exist, of course, and it appears likely that the 
variation of the absorption with time is, to a large extent, the result of 
varying C3 number density in the viscous region. Therefore, this effect is 
minimized by evaluating the absorption at the shortest possible time after 
shock passage. 

Accordingly, the present data were evaluated 50 psec after shock passage. 
The temperature and pressure at this time, designated T^q and p^Qi were 
2 percent and 5.5 percent, respectively, greater than the postshock values. 

The number densities of C3 used for the present data were obtained by first 
calculating the number density immediately behind the shock wave with the 
equilibrium chemistry program ACE and then using curve-fit relations as a per- 
turbation on this value to account for the change in composition as the pressure 
and temperature rose from the values P2 and T2 to P50 and Tgg* 

Spectrum Identification 

Since a number of gas species are expected to be present in the test-gas 
sample, the portion of the total absorption spectrum which is attributable to 
C3 must be determined. Figure 10 shows the mole fractions of the principal 
species which are predicted to be present behind the incident shock wave as a 
function of shock velocity. For this example, the initial pressure was held 
constant at 1.72 kPa, which is typical for the present data. Argon, which of 
course is the most abundant species, has been omitted from figure 10. The fig- 
ure shows other important species present to be H, H2, C, C2, C3, C2H, and C3H. 
Of these, the first four have well-identified spectra which may be recognized 
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and subtracted out if they appear in the present data. The spectra of the lat- 
ter two, C2H and C3H, are not known; therefore, any contribution they may make 
to the present absorption data could go unrecognized. However, both species 
have distinctly different variations with temperature than C3, and thus a com- 
parison of high-temperature and low-temperature absorption spectra should be 
helpful in identifying any contribution these species may make. In refer- 
ence 19, a spectrum for C2H was synthesized by analogy to similar molecules. 
This synthesized spectrum does not include any transitions in the wavelength 
range of the present data. 


RESULTS AND DISCUSSION 

The basic experimental results are presented in table II for the 3-percent 
acetylene mixture and table III for the 2-percent mixture. The shock velocity 
and postshock temperature and pressure T50 and P50 are tabulated as well as 
the absorption cross sections obtained from the eight spectrograph channels. 

For each channel, the central wavelength for each light pipe is noted, and the 
corresponding cross-section measurement is given directly below. For certain 
wavelength settings of the grating, some "stray light" was noted on channels 1 
and 2. When this occurred, data from these channels were not included in the 
tables. Here the term "stray light" refers to light emanating from the pulsed 
light source, but of the wrong wavelength for that particular channel. At each 
grating position, a check for stray light was made by means of a series of 
band-pass filters inserted in the light beam at the entrance slit to the 
spectrograph. 


Constant-Temperature Survey 

Figure 1 1 presents the absorption cross section for seven runs in which 
the postshock test state was repeated as closely as possible. The temperature 
was 3723 ± 47 K for these runs, which are runs 1 to 7 in table II. Between 
these runs, the grating position was changed to move the spectral setting of 
the channels by about 35 nm. Thus, a complete spectral survey was obtained at 
nearly constant test conditions. Figure 11 shows the photomultiplier data for 
the seven runs as square symbols. For two of these runs, the data from the 
densitometer readings of the film are shown as solid lines. The data from the 
film are seen to be in good agreement quantitatively with the photomultiplier 
data, particularly near the absorption peak, and confirm that the absorption 
profile is essentially smoothly varying, with the only distinctive features 
being due to the contaminant CN and the band heads of the C2 Swan system. 

In order to compare the present results with those of Brewer and Engelke 
(ref. 8), another survey was made at lower temperature. Figure 12 shows this 
survey, for which the mean temperature was 3240 K, and the results of refer- 
ence 8, for which the mean temperature was 3200 K (solid curve). To obtain 
this latter curve, the absorption data shown in figure 1 of reference 8 were 
used, and the number density of C3 was computed by using the thermal data of Lee 
and Sanborn (ref. 15) and AHg = 198 kcal/raol. This results in a computed C3 
number density of 8.89 x 10^5 molecules/cm^ for the Brewer and Engelke experi- 
ment, whereas a number density of 2.15 x 10^5 molecules/cm3 was obtained by 
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Brewer and Engelk^ using the free-energy functions of Pitzer and Clementi 
(ref. 20) and AHq = 188 kcal/mol. Figure 12 shows that use of the revised 
number density brings the furnace data at 3200 K of reference 8 into good agree- 
ment with the present shock-tube data at T = 3240 K. It should be noted that 
the violet band system of CN and the Swan band system of C£ were also observed 
by Brewer and Engelke but were "faired" out. (See fig. 1 of ref. 8.) 

The data at 3240 K represent the lower temperature limit for the present 
experiment, since the equilibrium density program predicts that solid carbon 
will be present in the flow for lower temperatures. No effects of particulate 
carbon are observable in the data of figure 12. A discussion of tests designed 
to observe particulate carbon is contained in a subsequent section of the 
report . 

No data were presented in reference 8 for wavelengths shorter than 370 nm. 
For less than 370 nm, the present data (figs. 11 and 12) show a continuation of 
the downward trend to about 300 nm, at which point the cross section becomes 
nearly constant (fig. 11). The data for this shortest wavelength region (below 
300 nm) must be viewed with some skepticism. The light output of the xenon 
flash lamp was very low in this region, resulting in poor signal-to-noise ratio 
on the photomultiplier channels and susceptibility to stray light contamination. 
Also, insufficient data were obtained at these wavelengths to indicate whether 
the absorption was due to Cg or some other species. 


Temperature Variation of Absorption Cross Section 

For each setting of the spectrograph grating, a number of runs were made 
with varying shock strengths, and thus, the temperature effect on absorption 
was determined at a number of wavelengths. A total of 26 different wavelengths 
was available with which to determine the temperature effect on absorption cross 
section. In figure 13, for six of these wavelengths the experimentally deter- 
mined absorption cross section is shown over the temperature range of the exper- 
iment. Also shown for each wavelength is a least-squares curve fit to the data 
of the form 

= aQ + a-|T + a2T^ + a^T^ (8) 

The curve-fit equations may then be used to generate an approximate spec- 
tral absorption profile for any specified temperature within the temperature 
range of the experiment. Such profiles are presented for four selected temper- 
atures in figure 14, where the value of equation (8) is shown at each wavelength 
by a symbol. The vertical bar shows the standard deviation of the data from the 
curve for that wavelength setting. Absence of a bar indicates insufficient data 
to determine standard deviation. 

The spectral absorption curve fits shown in figure 14 indicate that for the 
temperature range 3400 to 4000 K, the profiles remain quite similar. However, 
there appears to be some change with temperature in the general amplitude. The 
amplitude is highest at the low temperature extreme (3400 K) and decreases to a 
minimum at about 3800 K. Figure 13 shows that this trend is displayed at most 
wavelengths, with the minimum usually occurring near 3800 K. At temperatures of 
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4000 K and above, the general trend is for the absorption cross section to 
rise, but the very limited data in this temperature range do not permit a firm 
conclusion. 


It can be shown (e.g., ref. 21) that for a given atomic transition, the 
integral of the spectral absorption coefficient is proportional to the elec- 
tronic oscillator strength. Thus, 


I 


do) = 


ve‘ 


mc^ 


Nfc 


(9) 


in which w is wave number, fg is the electronic absorption oscillator 
strength, and N is the number density of the absorbing state. Now, as is 
done frequently, the concept of oscillator strength will be assumed to apply 
to an entire molecular transition, even though there are weaknesses in making 
this assumption. Since, for the present transition, the absorbing state is the 
ground state, N may be taken as the total number density of C3. Equation (9) 
may then be solved for fg to give 


fe = 


1.13 X 10 
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For the present results, the integral of the absorption cross section was 
approximated by using the curve fits (eq. (8)) for each wavelength to define o^j 
as a function of temperature, and the integral was determined by using the trap- 
ezoid rule. For this calculation, certain wavelength channels were believed 
to be influenced by the C2 Swan bands, and therefore, only 21 channels were 
included in the computation. The results of this computation are presented in 
figure 15. This figure shows that the oscillator strength varies over the tem- 
perature range of the experiments from about 0.062 at 3300 K to a minimum of 
0.036 at about 3900 K. The average oscillator strength is 0.04? over the tem- 
perature range 3300 to 4200 K. 


An estimate of the possible errors in the oscillator strength indicates 
that the greatest uncertainty lies in the calculation of C3 number density in 
the flow. An uncertainty of ±40 percent is assigned to this calculation, which 
is combined with an estimated maximum systematic error in the absorption mea- 
surement of ±20 percent to yield an overall uncertainty of 60 percent. 

Brewer and Engelke (ref. 8) calculate an oscillator strength fg of 
0.13 for their experiment. However, when their C3 number density is recalcu- 
lated according to the pj-esent recommendations (Lee and Sanborn's thermal prop- 
erties (ref. 15) and AHq = 198 kcal/mol), the oscillator strength becomes 
0.031. The difference in this value and the value obtained in the present 
experiment is caused by the more extensive wave-number range over which the 
present absorption data are integrated. 
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Particulate Formation 


One of the major concerns in the present experiments is whether solid par- 
ticulates were present in the test stream. For the present test conditions, 
the ACE computer program predicts that solid carbon will be present if the 
postshock temperature is below about 3220 K. For this reason, the present 
tests used only runs for which the postshock temperature was computed to be 
3300 K or greater, with the exception of the low-temperature survey presented 
in figure 1 2 . 

If thermochemical equilibrium was not achieved in the present tests, and 
as a result, the test stream contained particulate carbon, the absorption 
results would include the absorption and scattering of the particulates. Also, 
the actual concentration of carbon-bearing gaseous species would be lowered 
because of the loss of carbon to the solid state. 


In order to examine the effects of particulate carbon in the test gas, a 
series of runs was made at postshock temperatures below 3220 K. Figure 16 
shows examples of the absorption coefficients measured. As the temperature of 
the postshock test gas is decreased, the character of the spectral absorption 
profile changes: the absorption coefficient increases markedly for the short- 

wavelength portion of the spectrum. Thus, the absorption coefficient rises 
monotonically with decreasing wavelength at temperatures below 3000 K 
(fig. 16(c)); whereas at temperatures above 3200 K (e.g., fig. 16(a)), a max- 
imum occurs near 400 nm, with decreasing absorption at shorter wavelengths. 

The calculated effects of particulate absorption and scattering are also 
shown in figure 16. The absorption and scattering were calculated by the 
method described in reference 22. According to reference 22, most soot parti- 
cles are either "spheres" or long strands with radii lying between 10 and 
40 nm, and thus 2iTa/X « 1, where a is the particle radius. For this 
"Rayleigh region" the absorption cross section c^X^abs given as 

‘^X,abs = 

where pj^ is the so-called "Rayleigh factor." The absorption coefficient is 
then 

3 m 

^X,abs = Np'^X,abs " "h 7” ^X 

4 Pc 

where m is the mass of particulates per unit flow volume and is the 

density of carbon. Interpolated values for pj^ as a function of wavelength 
were obtained from the tabulated values given in reference 22. 

The Rayleigh scattering cross section is given in reference 11 as 


2 

24ir3v^ - N2^ - 1)(Ni2 - N2^ + 2) + 4 Nt2n22] + 36Ni2n22 

X^ [(Ni2 _ N 22 + 2)2 + 4N-|2n22]^ 
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where V is the particle volume and N-| + IN2 is the complex refractive index 
of the carbon particles. Values for N-|^ - ^ 2 ^ and 2N-]N2 were obtained from 
reference 23. The contribution to the extinction coefficient due to Rayleigh 
scattering is then given by 


3 m 

^XjSca = ^p*^X ,sca "1. _ *^X ,sca 

^ PcTra^ 

In figure 16, the total contribution to the extinction coefficient due to par- 
ticle absorption and scattering is shown for an assumed particle radius a of 
30 nm. Since the absorption term dominates and is independent of particle 
radius, the resulting total is not very sensitive to assumed particle size. 

The calculated extinction coefficients have been reduced by a factor of 4 in 
order to compare them better with the experimental data. It is seen in fig- 
ure 16 that the change in the nature of the experimental absorption profile as 
the temperature is reduced below 3000 K may be ascribed to the effects of par- 
ticulate carbon. The fact that the extinction coefficient is lower by about a 
factor of 4 than the calculated level suggests that less particulate carbon is 
present in the experiment than is calculated. These results strongly suggest 
that no appreciable amount of particulate carbon is present at the higher tem- 
peratures for which the present experimental program was conducted. 


CONCLUDING REMARKS 

The spectral absorption properties of triatomic carbon C3 have been mea- 
sured in a shock tube containing a test-gas mixture of acetylene diluted with 
argon. The temperature behind the incident shock ranged from about 3240 to 
4300 K, and the corresponding pressure varied from 37 to 229 kPa. Appreciable 
absorption was measured at wavelengths from about 300 nm to about 540 nm. 
Absorption was also noted below 300 nm, but the lack of sufficient data and the 
poor signal-to-noise ratio in this wavelength range precluded a firm determina- 
tion that this short-wavelength absorption was due to C3. 

The absorption profile was very similar to that measured by Brewer and 
Engelke (J. Chem. Phys., vol. 36, no. 4, Feb. 1962) in a furnace at 3200 K. 

They did not make measurements below 370 nm and thus did not include the short- 
wavelength portion of the profile determined in the present study. However, 
for the overlapping range of wavelengths of the two experiments, the results 
are in good agreement . 

Data taken for a narrow passband at specific wavelengths over a range of 
temperature and pressure reveal a temperature variation of the spectral absorp- 
tion cross section. Curve fits for these data were obtained for 21 wavelengths, 
and these curve fits can be used to construct approximate spectral profiles for 
any selected temperature within the valid range of the data. If such profiles 
are integrated with respect to wave number, an electronic oscillator strength 
is obtained for the transition. This oscillator strength varied from 0.062 to 
0.036 over the temperature range 3300 to 4200 K and may be compared with a value 
of 0.13 quoted by Brewer and Engelke or a value of 0.031 if the number density 
of their experiment is reevaluated. 
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The absorption cross section is determined not only by the absorption 
coefficient measured, but also by the computed number density of the absorbing 
state. This latter quantity was computed in the present study by use of an 
equilibrium, free-energy-minimization computer program which assumed a value 
of the heat of formation for C3 of 198 kcal/mol at a temperature of 0 K. The 
thermal properties of C3 were taken from the calculations of Lee and Sanborn 
(High Temp. Sci., vol. 5, no. 6, Dec. 1973). 

Solid carbon particulates were not expected to be present at the postshock 
temperatures of the present tests, on the basis of the equilibrium calculations. 
In order to examine the effects of particulates, a few runs were made at lower 
temperatures, where solid carbon was calculated to be a part of an equilibrium 
composition. As the temperature was decreased below 3000 K, the absorption 
profile underwent a distinct change and assumed a shape similar to that which 
is calculated for absorbing and scattering Rayleigh particles. The experimental 
profile was, however, lower than the calculated profile by about a factor of 4, 
suggesting that less carbon was in the solid state than was calculated. It was 
concluded, therefore, that in the present tests, the test gas at postshock tem- 
perature above 3300 K was free of solid carbon particulates. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 1, 1978 
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THERMAL PROPERTIES OF C 3 

The determination of the concentration of C 3 behind the shock wave for 
the present test conditions assumes that thermochemical equilibrium has been 
attained in the flow, and consists of minimizing the Gibbs free energy for the 
mixture of species which are present. This calculation requires a definition 
of the temperature-dependent thermodynamic properties for each species. In 
order to have a common reference point for all species, the enthalpy is 
referred to the elements in their most natural form at 298 K and a pressure 
of 1 atm. In this form, the equation for static enthalpy is 

pT 

h= Cp dT + AH 298 (A1) 

^298 

O 

Equation (A1) shows that the heat of formation AH 298 is required for 
each species vrtiich is to be considered in the final mixture. In this appendix 
the available information on the thermodynamic properties and heat of formation 
of C 3 is reviewed in order to assess the uncertainty in calculated number den- 
sities. A thorough review of the state of knowledge of carbon vaporization as 
of 1969 is given in reference 24. This appendix attempts a much briefer review, 
but includes data and calculations published since 1969. 


Evaluation of Vapor-Pressure Data 

The heat of formation of C 3 has not been established to a satisfactory 
degree of accuracy, despite the number of experiments which have b- 3 en conducted 
to determine this quantity. The experiments are based on the vaporization of 
graphite. The chief experimental difficulty is the determination of the frac- 
tion of total vapor pressure which is due to the partial pressure of C 3 . 

The vapor-pressure data are commonly einalyzed by two methods, referred to 
as the second-law method and the third-law method. The second-law method makes 
use of the Clausius-Clapeyron equation, with assumptions that the vapor behaves 
as an ideal gas and the changes in specific volume of the solid are negligible 
(see, for example, ref. 25): 


AHf 


d(ln p) 

— R 

d(1/T) 


(A2) 


The logarithm of the experimental values of the partial pressure of C 3 is' plot- 
ted against the inverse of the vapor temperature, and the heat of formation may 
thus be inferred from the slope of the line through the data. 

The advantage of the second-law method lies in the fact that only relative 
values of pressure, not absolute values, are required. Since the partial pres- 
sure of C 3 is usually obtained from mass spectrometer readings, the scaling 
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factors which are used to convert the spectrometer signal of the mass peak 
to partial pressure do not influence the slope of the curve. The chief disad- 
vantage of this method is that the equation yields the heat of formation at the 
temperature T of the pressure measurements. Now if only a small range of 
temperature is used in the experiments, an accurate determination of the slope 
of the line through the data is difficult. If a wide range of temperature is 
used in the experiment, the data may not be expected to be along a straight 
line, since the heat of formation is, in general, temperature dependent. 

The heat of formation at the temperature T may be converted to the value 
at the reference temperature by Kirchhoff's equation: 

O O 

AHj,ef = AHx + (hx - hi,gf)goiid - O^x - hr.ef)yapor (^3) 

where the subscripts ref and T denote the reference temperature and the 
given temperature. Both 0 K and 298.15 K have been used as the reference tem- 
perature. However, data based on one reference temperature may be converted to 
another reference by equation (A3). Note that equation (A3) requires a knowl- 
edge of the enthalpy-temperature variation for both the vapor and solid states. 

The third-law method makes use of the free-energy function and may be 
written 


AHref 

T 


-R In p + 


Gt - hr^ef) 


solid 


Gt - hr^ef ' 


vapor 


(A4) 


Reference 26 discusses the evaluation of vaporization data by the second-law 
and third-law methods. Although the third law allows a determination of the 
heat of formation from a single measurement, it requires a knowledge of the 
Gibbs free energy for both the vapor and solid states. Furthermore, the cal- 
culation now depends upon an absolute determination of the species vapor 
pressure p. 


It is common practice to evaluate a set of data by both the second-law and 
the third-law methods, to note the consistency of results. However, the results 
of the third-law method are to be preferred over those of the second-law method, 
because evaluation of each data point is independent of the others. This allows 
trends with temperature or inconsistent data points to be noted more easily. 

The foregoing discussion has reviewed these data evaluation methods in 
order to emphasize the importance of an accurate knowledge of the thermodynamic 
properties of the species under study. Specifically, the enthalpy and Gibbs 
free energy must be known over the temperature range of interest, and if a 
change is made in the assumed variation of enthalpy and Gibbs energy with tem- 
perature, the data will yield a different value of heat of formation. 


Thermodynamic Functions of C3 

The free-energy function for was calculated by Glockler (ref. 27) in 
1954. He assumed a multiplicity of 5 for the ground state and a bending f re- 
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quency of 353 cm“^ and estimated the heat of formation as 212.3 kcal/mol. His 
calculation was based on the standard rigid-rotor harmonic-oscillator model. 

In 1959 Pitzer and Clementi (ref. 20) corrected the multiplicity of the ground 
state to 1 and recomputed the free-energy function and the enthalpy using a 
bending frequency of 550 cm“^. However, in 1965 Gausset , Herzberg, Lagerqvist, 
and Rosen (ref. 3) published a detailed spectroscopic analysis of C 3 , conclud- 
ing that the bending frequency of the ground electronic state was of the order 
of 65 cm“^ for the first few levels. This low bending frequency created a 
difficulty in calculating the thermal functions , for as shown by Robiette and 
Strauss (ref. 28), the assumption of a harmonic oscillator with the lower value 
of bending frequency causes the calculated entropy to be too high compared with 
that obtained from the temperature measurements. For example. Thorn and Winslow 
(ref. 29 ) calculated a value of entropy S 2400 77.^1 eu from their vapor- 

pressure measurement, and Robiette and Strauss (ref. 28) interpreted the pres- 
sure measurements of Drowart et al. (ref. 13) to yield S 2400 = 76.1 eu. 

Robiette and Strauss showed, however, that at 2400 K a harmonic-oscillator 
calculation with a bending frequency of 63 cm~^ yields an entropy of 85.56 eu. 
They showed that the entropy can be reduced by about 3 eu by taking the quartic 
anharmonicity into account . 

Strauss and Thiele (ref. 30) calculated the thermal properties of C 3 using 
the spectroscopic data of Gausset et al. (ref.. 3) and accounting for rotation- 
vibration interaction. However, as they pointed out, only the first few levels 
of the bending vibration have been measured, so that the potential function 
must be based on a very limited set of energy levels. They show that inclusion 
of the quartic anharmonicity lowers the entropy at 2400 K by about 2.8 eu, and 
inclusion of the vibration-rotation coupling lowers it by another 1.1 eu, to a 
value of 81 . 4 . Further, they demonstrated that even if the potential increased 
drastically above a bending angle of 0.4 tt, the calculated entropy cannot be 
lowered below S 2400 = 79.8 eu. It appears therefore, that the best calcula- 
tion method cannot compute entropy values as low as those inferred from the 
experiments of references 13 and 29 . 

The JANAF tables (ref. 26) adopted a potential function which makes use 
of the vibrational frequencies measured by Gausset et al. (ref. 3) for the low 
levels and then steepens to fit a vibrational frequency interval of 650 cm"^ 
for higher levels. The value of 650 cm~^ is admitted to be arbitrary and was 
selected to yield functions which approximate experimental data. Their value 
of entropy at 2400 K is S 2400 “ 78.4 eu. 

Hansen and Pearson (ref. 31) used a quadratically perturbed square well 
potential in a quantum calculation to approximate the spectroscopic data of C 3 . 
The maximum molecular bending at the potential wall was 64° for this model, 
which is in close agreement with Strauss and Thiele's approximate model 
(ref. 30 ), which steepened the potential at a bending angle of 0.4 tt rad, or 
72 °. From their analysis, Hansen and Pearson calculated an entropy at 2400 K 
of S 2400 = 79.68 eu, and concluded that their limiting calculation thereby 
supports the calculation of Strauss and Thiele. 

Lee and Sanborn (ref. 15) also calculated the thermal properties of C 3 by 
developing a special potential model. Their enharmonic model for the lower 
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vibrational levels fits the data of Gausset et al. (ref. 3), and for higher 
vibrational levels the term difference approaches 308 cm“^, the bending 
frequency of the first excited state. They also included the electronic 
contribution to the partition function, which had been neglected in the other 
treatments. However, they used a rigid-rotor assumption and thus neglected 
vibration-rotation interaction. Lee and Sanborn showed that their results are 
in close agreement with those of Strauss and Thiele (ref. 30) when they omitted 
the electronic contribution. They used the electronic energy levels given by 
Pitzer and Clementi (ref. 20) in computing the electronic contribution to the 
partition function. 

Figure 17 compares the specific heat at constant pressure, the enthalpy, 
the entropy, and the Gibbs energy function as calculated by references 15, 26, 
30, and 31. The differences between the calculations are most accentuated in 
Cp, but are appreciable in the other functions as well. The Gibbs functions 
as computed by Lee and Sanborn are seen to be quite close to those by Strauss 
and Thiele throughout the temperature range. About half the difference between 
these curves may be attributed to the electronic contribution which Lee and 
Sanborn included. 

The calculations of Strauss and Thiele (ref. 30) would appear to be the 
most rigorous treatment of the C3 molecule with the spectroscopic data avail- 
able, except for omission of the electronic contribution to the partition func- 
tion. Since the results of Lee and Sanborn are in acceptable agreement with 
the results of Strauss and Thiele when the electronic contribution is omitted, 
their thermal functions (ref. 15) seem to be the best choice. Lee and Sanborn 
also presented thermal functions for C4 through Cy. In addition, the proper- 
ties for C-| and C2 were extended to 10 000 K. Thus a consistent set of prop- 
erties for the carbon species up to Cj are available in reference 15 for carbon 
vapor calculations. 

For the reasons stated, the equilibrium calculations of this study used 
curve fits to the calculated data of Lee and Sanborn for the thermodynamic 
functions of C3. 

One further point of interest may be noted regarding the evaluation of 
the various calculations of thermal properties. The most nagging problem has 
been the fact that each of the calculations yielded higher entropy values than 
those measured by Drowart et al. (ref. 13) and Thorn and Winslow (ref. 29), and 
resolving the low bending frequency with the experimental entropy values seemed 
impossible. In 1972, however, Wachi and Gilmartin (ref. 32) made new vapor- 
pressure measurements with a calibrated mass spectrometer and a Knudsen effu- 
sion cell. They quoted experimental entropies at 2400 and 2700 K as 
S2400 = 81.47 eu and S2700 = 82.91 eu. These values are shown in fig- 
ure 17(c) and are in excellent agreement with the calculation of Strauss and 
Thiele (ref. 30) and in acceptable agreement with the calculation of Lee and 
Sanborn (ref. 15). Thus, if one assumes that the measurements of Wachi and 
Gilmartin are superior to the earlier ones, the discrepancy between theory and 
experiment is resolved on this point. 
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Heat of Formation of C 3 

In the previous section, reasons were given for choosing the calculated 
data of Lee and Sanborn (ref. 15) for curve fitting the equilibrium calcula- 
tions of the present study. Having made a choice, the various experimental 
vapor -pres sure measurements may now be compared using the same thermal property 
data for each. In addition to the measurements of Drowart et al. (ref. 13), 
mass spectrometer measurements have also been made by Wachi and Gilmartin 
(ref. 32 ), by Milne et al. (ref. 33), and by Zavitsanos and Carlson (ref. 34). 

The conversion of mass spectrometer signal level to partial pressure 
requires several conversion factors relating to ionization cross section, sec- 
ondary electron yields of the electron multiplier, and spectrometer transmis- 
sion coefficient. In reference 35, Meyer and Lynch made a careful evaluation 
of the factors affecting spectrometer response, recommended new conversion fac- 
tors, and reevaluated the data of references 13 and 34. They also included the 
data of reference 32, but since Wachi and Gilmartin used gold and cobalt to 
calibrate their mass spectrometer, their conversion factors were not revised. 
Meyer and Lynch then used the thermal properties of Strauss and Thiele 
(ref. 30 ) to reduce the data. 

In the present study the conversion factors of Meyer and Lynch (ref. 35) 
were adopted and applied to the data of references 13, 33, and 34. The data of 
Wachi and Gilmartin (ref. 32) were used without correction. The resulting par- 
tial pressure of C 3 is shown in figure I 8 as a function of the inverse tempera- 
ture 1/T. 

It was deemed highly desirable to include data obtained by a different 
method, and in a different temperature range; therefore, the data of Lundell 
and Dickey (ref. 36) have also been included in figure I 8 . Lundell and Dickey 
used a powerful laser to vaporize graphite, and inferred the total vapor pres- 
sure from the surface recession rate. Figure I 8 includes the five data points 
for which Lundell and Dickey state that the condition of pure carbon vapor at 
the surface was achieved. The partial pressure of C 3 was computed by subtract- 
ing the calculated partial pressures of C and C 2 from the total vapor pressure. 
The contributions of Cij and larger molecules are expected to be negligible; how- 
ever, note that neglect of these species means that the pressure shown in fig- 
ure 18 for these data is an upper limit for the partial pressure of C 3 . 

All the vapor-pressure data shown in figure I 8 were reduced to heat of 
formation by using the third-law method (eq. (A4)) and the thermal properties 
of Lee and Sanborn (ref. 15). For the thermal properties of solid state car- 
bon, the curve fits of Gordon and McBride (ref. 37), which fit the values given 
in the JANAF tables (ref. 26), were used. 

^ The average value of heat of formation for all data shown in figure I 8 is 
AHq = 197.85 kcal/mol. However, it may be argued that this method of averaging 
tends to bias toward the experiments with the most data points; therefore, the 
data were also averaged by experimenter. The^^ results are shown in the table on 
the following page and give an average of AHq = 199.50 kcal/mol. Determina- 
tion of a second-law value does not appear to be helpful. Since the heat of 
formation varies with temperature, and since the composite data cover a wide 
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Investigators 

Average AHq , 
kcal/mol 

Drowart, Burns, DeMaria, and Inghram (ref. 13) 

Wachi and Gilmartin (ref. 32) 

Zavitsanos and Carlson (ref. 3^) 

Milne, Beachey, and Greene (ref. 33) 

Lundell and Dickey (ref. 36) 

Average 

Overall average of all data points 

197.25 

200.84 

198.73 

196.22 

204.48 

199.50 

197.85 


temperature range, a single straight-line fit is not meaningful. Moreover, the 
local slope of a curve fit through the. data is very dependent on the analytical 
form of the assumed curve. In order to illustrate the curvature of a constant 

O 

AHq cyrve, three computed curves are shovm in figure 18 for constant values 
of AHq = 19^, 198, and 202 kcal/mol. 

The conclusion drawn from this reexamination of carbon data is that the 
range of values for the heat of formation of which may be defended is from 
197.85 to 199.50 kcal/mol. A value of 198 kcal/mol was selected for the pres- 
ent study and Incorporated into the computer code which was used to determine 
C3 number densities for the test conditions. 


2H 
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TABLE I.- SHOCK-TUBE INSTRUMENTATION PORT LOCATIONS AND TYPE OF DETECTOR 


Station 


11 (test window) 


^Abbreviation for type of detector: 
H.T. thin-film heat-transfer gage 
P. pressure transducer 
P.M, photomultiplier 


Distance from 
diaphragm, m 

Detector 

(a) 

2.889 

H.T. 

4.417 

P. 

5.925 

P. 

7.447 

H.T. 

8.357 

P.M. 

8.662 

P. 

9.566 

P. 

11.092 

P. 

12.617 

H.T. 

14.126 

P. 

14.940 

P.M. 

15.149 

P. 










TABLE II.- TEST CONDITIONS AND MEASURED ABSORPTION CROSS SECTIONS 
FOR TEST -GAS MIXTURE: 2.99 PERCENT C 2 H 2 , 97.01 PERCENT ARGON 



If 

For each spectrograph channel and each run, X is given on the top line; 
a, on the bottom line. 
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TABLE II.- Continued 


^s» T50, P 50 , (NCo)cq X 10 

Run m/sec K kPa ' ' 


X, nm, and CJ;\^ x 1 q 1®, cm^,* 
for spectrograph channel - 



For each spectrograph channel and each run, X is given on the top line; 
a, on the bottom line. 
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TABLE II.- Continued 


Run 

Vs. 

m/sec 

T50. 

K 

P50. 

kPa 

(“=3)50 X 

X, nm, and x 10^®, cm^,* 

for spectrograph channel - 

1 

2 

3 

4 

5 

6 

7 

8 

21 

2317 

3896 

128 

2.60 



405 

434 

461 

490 

520 

545 



— 

— 

3.68 

1 .55 

0.98 

0.67 

1.57 

1.13 


2429 

4103 

102 

1 .58 



405 

434 

461 

490 

520 

545 



— 

— 

3.64 

2.05 

2.24 

— 

2.69 

2. Oil 

23 

2288 

3721 

37 

0.86 



405 

434 

461 

490 

520 

545 



— 

— 

2.95 

1.19 

0.82 

0.65 

0.26 

0.49 

24 

2117 

3411 

53 

1 .06 



405 

434 

461 

490 

520 

545 





4.97 

1 .82 


0.61 

0.27 

0.57 




25 

2096 

3380 

62 

1.12 



405 

434 

461 

490 

520 

545 





4.49 

1 .71 

0.49 

0.40 

0.29 

0.43 



26 

2573 

4328 

65 

0.43 



405 

434 

461 

490 

520 

545 




— 

7.43 

6.23 

5.45 

2.71 

4.21 

5.99 

27 

2506 

4224 

79 

0.85 



405 

434 

461 

490 

520 

545 






3.29 

2.96 

1 .47 

1 .06 

3.50 




28 

2434 

4112 

99 

1.52 



405 

434 

461 

,490 

520 

545 





4.83 

2.47 

2.06 

0.91 

1.65 

2.23 



29 

2237 

3721 

119 

2.48 



405 

434 

461 

490 

520 

545 





4.66 

1 .70 

0.98 

0.76 

0.46 

1.10 



30 

2243 

3759 

178 

3.56 



405 

434 

461 

490 

520 

545 



— 

— 

4.09 

1.66 

0.85 

0.37 

0.45 

0.78 


*For each spectrograph channel and each run, X is given on the top line; 
a , on the bottom line. 
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TABLE II.- Continued 



*For each spectrograph channel and each run, X is given on the top line; 
a , on the bottom line . 
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TABLE II.- Concluded 


■ 


Run 

— 

Vs, 

m/sec 

T50, 

K 

P50, 

kPa 

(Nc3)5o X 10-16 

X, nm, and a \ x 10^®, cm^,* 
for spectrograph channel - 

1 

2 

3 

4 

5 

6 

7 

8 

41 

2221 

3666 

116 

2.37 

— 

— 

256 

272 

285 

299 

314 

328 

— 

— 

2.60 

2.00 

0.90 

1 .40 

1 .50 

1 .78 

42 

2309 

.3845 

89 

1.86 

— 

— 

256 

272 

285 

299 

314 

328 

— 

— 

1 .70 

1 .58 

1 .24 

0.80 

1 .54 

1.56 

43 

2082 

3370 

59 

1 .01 

— 

— 

256 

272 

285 

299 

314 

328 

— 

— 

3.20 

2.50 

1 .60 

1.51 

2.46 

2.30 

44 

2038 

3264 

59 

0.864 

287 

299 

314 

329 

343 

357 

371 

385 

— 

— 





6.37 

6.46 

45 

2035 

3260 

58 

0.854 

325 

338 

352 

368 

• 

381 

395 

409 

423 

2.67 

4.10 

5.10 

6.03 

6.04 

6.62 

4.33 

2.35 

46 

2026 

3242 

58 

0.807 

360 

373 

388 

403 

416 

430 

445 

458 

6.12 

6.30 

6.97 

6.21 

— 

2.35 

1.03 

0.58 

47 

2026 

3242 

59 

0.823 

431 

444 

458 

474 

487 

501 

515 

529 

2.06 

0.75 

0.76 

1.60 

0.30 

0.40 

2.10 

0.13 

48 

2011 

3212 

57 

0.748 

395 

408 

423 

438 

452 

465 

480 

494 

6.88 

4.45 

2.39 

1.73 

0.70 

1.58 

0.57 

0.56 

49 

2017 

3221 

57 

0.769 

466 

479 

494 

509 

522 

536 

551 

564 

1 .87 

0.79 

0.52 

1 .71 

0.64 

0.63 

0.67 

1 .40 


For each spectrograph channel and each run, X is given on the top line; 
a , on the bottom line. 
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TABLE III.- TEST CONDITIONS AND MEASURED ABSORPTION CROSS SECTIONS 


FOR TEST -GAS MIXTURE; 1.97 PERCENT C 2 H 2 , 98.03 PERCENT ARGON 


T50, P50, (NCojcQ ^ 
Run m/sec K kPa ^ 



X , nn 
for 

1 , and x 10^®, cm^,* 

spectrograph channel - 


1 


3 

D 

5 

6 

m 

8 

B 




Bji 


2.00 

0.93 

1 .75 

0.30 

m 

2.30 

BSS' 

my 


B 

520 

545 


— 

2.00 

0.73 

1.20 

— 

1.39 

B 

— 

405 

434 

461 

490 

Qymni 

B 



0.30 



57 

2052 

3525 

195 

2.25 



0.24 

0.10 

0.94 

490 

520 

545 

0.27 

0.32 

0.70 




0.14 

0.31 

0.49 


1 

490 

520 

545 

0.60 

0.25 

00 

00 

0 

0.63 

461 

490 

520 

545 

0.88 

0.71 

0.47 

0.55 

461 

490 

520 


0.77 

0.34 

0 

0.56 


For each spectrograph channel and each run, X is given on the top line; 
a, on the bottom line. 













































































Xenon light source 


Quartz lenses- 


1.8-mm-diameter orifice 


-Quartz windows 


Shock tube 


Flow 

- 1 -cm -diameter field stops 


Dump tank 


-Mask 


Mirror- 


Light pipes 

and photomultipliers- 


-Spectrograph 


Mirror partially^ 
into beam 


■Drum camera 


Figure 1.- Schematic of shock-tube test location and optical arrangement. 







Increasing 

light 


Shock Interface 

arrival arrival 


(a) Photomultiplier traces. 


137 kPa 

I 


Station 3 
Station 4 



200 Msec 


(b) Sample pressure transducer traces. (Note that the traces 
have different start times.) 


L-78-15 


Figure 3.- Sample oscillograms showing photomultiplier and pressure transducers 
as shock detectors. P 2 = 221 kPa; T 2 = 3638 K. 
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Photomultiplier signal. 



t, ijsec 


(c) Digitized records of run and prerun traces. 



t, psec 


(d) Absorption parameter F(t). 

Figure 5.- Sample spectrograph channel record and method of determining 
absorption parameter F(t). pcQ = 105 kPa; Teg = 3709 K. 




Shock 

arrival 


Light pulse 
start 


a> 

■a 



Wavelength, nm 


l-78-17 


(a) Print of drum camera film 


Before shock 
arrival 


100 Msec 

after shock arrival 


320 340 360 380 400 420 440 460 

Wavelength, nm 


40 


(b) Microdensitoraeter tracings. 

Figure 6.- Sample drum camera record and microdensitometer tracings 

P50 = 118 kPa; T50 = 3686 K. 





Density above gross fog 
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,d shock 



(a) 

P 2 = 209 kPa; 

(b) P 2 = 99 kPa; 

(c) P 2 = 

94 kPa 


T2 = 3562 K; 

T2 = 3639 K; 

T2 = 

4026 K 


X = 43^ nm. 

X = 317 nm. 

X = 

461 nm. 

Figure 9-- 

Example of time variation 
after correction 

of absorption F(t) and resulting 
for free-stream state variation. 

cross section 

c^X 




Figure 10.- Mole fractions of principal species behind a shock wave into 
a gas mixture of 97 percent argon, 3 percent C 2 H 2 . Pi = 1.72 kPa. 

( Argon is not shown . ) 


















(a) T = 3400 K. 



(b) T = 3600 K. 

Figure 14.- Spectral cross-section profile at selected temperatures 
Symbols are value of cubic curve fit; bars indicate standard 
deviation of data (where sufficient data available). 





r> 








Wavelength, X,nm 


Figure 16.- Measured extinction coefficient at various temperatures 
compared with calculated extinction due to scattering and 
absorption of 30-nm-radius particles. 



(a) Specific heat at constant pressure. 

Figure 17.- Comparison of thermal properties of C 3 
as calculated by four methods. 
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(d) Gibbs function. 
Figure 17.- Concluded. 
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